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Abstract. For the very first time, the correlation femtoscopy method is applied 
to a kincmatically complete measurement of meson production in the collisions 
of hadrons. A two-proton correlation function was derived from the data for 
the pp —> ppX reaction, measured near the threshold of n meson production. 
A technique developed for the purpose of this analysis permitted to establish 
the correlation function separately for the production of the pp + n and of the 
pp + pions systems. The shape of the two-proton correlation function for the ppn 
differs from that for the pp(pions) and both do not show a peak structure opposite 
to results determined for inclusive measurements of heavy ion collisions. 
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1. Introduction 



Momentum correlations of particles at small relative velocities are widely used to study 
the spatio-temporal characteristics of the production processes in relativistic heavy 
ion collisions pQ. This technique, called after Lednicky correlation femtoscopy [2J, 
originates from photon intensity interferometry initiated by Hanbury Brown and 
Twiss [3 . Implemented to nuclear physics [2j|4|[5] it permits to determine the duration 
of the emission process and the size of the source from which the particles are 
emitted [2J. A central role plays the correlation function which has been defined as 
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the measured two-particle distribution normalized to a reference spectrum obtained 
by mixing particles from different events [2J. The importance of the correlation 
femtoscopy has been well established for investigations of the dynamics of heavy ion 
collisions with high multiplicity. However, as pointed out by Chajecki [6], in the case 
of low-multiplicity collisions the interpretation of correlation function measurements is 
still not fully satisfactory, especially in view of the surprising observation by the STAR 
collaboration indicating universality of the resulting femtoscopic radii for both, the 
hadronic (proton-proton), and heavy ion collisions [7]. One of the challenging issues in 
this context is the understanding of contributions from non-fcmtoscopic correlations 
which may be induced by the decays of resonances, global conservation laws [6], or 
by other unaccounted interactions. In contrast to heavy ion collisions, in the case 
of single meson production, the kinematics of all ejectiles may be entirely determined 
and hence a kinematically complete measurement of meson production in the collisions 
of hadrons gives access to complementary information which could shed light on the 
interpretation of the two-proton correlations observed in heavy ion reactions. It is 
also important to underline that the correlation of protons was never exploited till 
now in near threshold meson productions, and as an observable different from the 
distributions of cross sections, it may deepen our understanding of the dynamics of 
meson production. Particularly favourable are exclusive experiments conducted close 
to the kinematical threshold where the fraction of the available phase-space associated 
with low relative momenta between the ejectiles is large [8]. 

In this article we report on a rj meson and multi-pion production experiment 
in which the mesons were created in collisions of protons at a beam momentum of 
2.0259 GeV/c corresponding to an excess energy of Q = 15.5 MeV for the pp —> ppr\ 
reaction. The measurement of the two-proton correlation function for these reactions 
is important not only in the context of studying the dynamics underlying the heavy 
ion physics. Such investigations are interesting by themselves because they offer a 
new promising diagnostic tool, still not exploited, for studying the dynamics of meson 
production in hadron collisions. 

The correlation function carries information about the emitting source and, in 
particular, about the size of the interaction volume of the pp — »■ pprj process. The 
knowledge of this size might be essential to answer the intriguing question whether the 
three-body pprj system is capable of supporting an unstable Borromean bound state. 
Borromean systems may be realized in a variety of objects on the macroscopic (e.g. 
strips of papers), molecular [9|fT0] and nuclear scale (e.g. n Li or 6 He nuclei [TTHT3] ). 
According to Wycech [TJ], the large enhancement of the excitation function for the 
pp —> pprj reaction observed close to the kinematical threshold may be explained by 
assuming that the proton-proton pair is emitted from a large (Borromean like) object 
whose radius is about 4 fm. 

2. Experimental technique 

The experiment was conducted using the proton beam of the cooler synchrotron 
COSY [15] and an internal hydrogen cluster target [16]. Momentum vectors of 
outgoing protons from the pp — > ppX reaction were measured by means of the COSY- 
11 facility [17] presented schematically in Figured] 

The two-proton correlation function R(k) was determined for the pprj and 
pp(pions) systems, respectively. Here, R(k) denotes a projection of the correlation 
function onto the momentum of one of the protons in the proton-proton center-of- 
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Figure 1. Schematic view of the COSY-11 detection setup |17| . Superimposed 
lines indicate trajectories of protons from the pp — ¥ ppr] reaction. The positively 
charged particles were identified using the time of flight between the scintillator 
detectors SI and S3 combined with a momentum reconstruction by tracking back 
the trajectories (reconstructed based on the signals from the drift chambers Dl 
and D2) through the magnetic field of the COSY dipole to the reaction point. The 
size of the detectors and their relative distances are not to scale. More detailed 
description can be found in | 17H20| . 

mass systenjjj. It was calculated, by means of the well established technique, as a 
ratio of the momentum (k) dependent reaction yield Y{k) to the uncorrelated yield 
Y*(k) according to the formula (cf. |21j ) 

where C denotes an appropriate normalization constant. Y*(k) was derived from the 
uncorrelated reference sample obtained by using the event mixing technique introduced 
by Kopylov and Podgoretsky [5]. 

The separation of the correlations of the pp + n from the pp + pions system and 
corrections for the limited acceptance of the detection system constitute the two main 
challenges to be solved when deriving the correlation function from the experimental 
data. 



2.1. Separation of events from the production of ppr\ and pp + pions systems 

In the discussed experiment, the four-momenta of the two final state protons were 
measured only, and the unobserved meson was identified via the missing mass 
technique [TH [20]. In such situations the entirely accessible information about the 
reaction is contained in the momentum vectors of the registered protons. Therefore, 
it is impossible to know whether in a given event the rj meson or a few pions have 

| Note, that some authors instead of k take as the independent variable the relative momentum of 
emitted particles q = \pi — pb.\ with k « q/2. 
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been created. However, statistically one can separate these groups of events using the 
missing mass spectra for each chosen region of the phase-space. 

Thus, Y pprj {k) can be extracted for each studied interval of k by grouping the 
sample of measured events according to the value of fc, next calculating the missing 
mass spectra of the pp —¥ ppX reaction for each sub-sample separately, and counting 
the number of pprj events from these spectra. Examples of typical missing mass 
histograms for three different k intervals are presented in Figure [2] The statistics 
obtained in the considered measurement allowed to divide the kinematically available 
range of k into bins whose width (Afc = 2.5 MeV/c) corresponds approximately to the 
accuracy of the determination of k with (a(k) « 2 MeV/c). A clear signal of ij meson 
production is observed in each spectrum on top of a continuous distribution originating 
from the multi-pion production. It is important to note that in the studied range of 
the missing mass the observed shape of the multi-pion distributions is well reproduced 
in corresponding simulation studies for each region of the phase space |18j . An 




Figure 2. Examples of missing mass spectra measured 1201 for the pp —¥ 
ppX reaction at k G (27.5; 30.0) MeV/c, k G (55.0; 57.5) MeV/c and 
k G (82.5; 85.0) MeV/c. Points represent experimental data. Dotted histograms 
depict simulations of the pp — » pp2ir, pp — > pp3ir and pp — ¥ pp4ir reactions, while 
solid line histograms present the sum of simulations for the pp — ¥ pprj reaction 
and for the pp —¥ pp2ir, pp —¥ pp3ir and pp — ¥ ppAir reactions, respectively. The 
simulated spectra were fitted to the data adjusting only the amplitude [20] . 

extraction of Y* prj {k) - unbiased by the multi-pion production - is, however, not trivial. 
Applying a mixing technique one can construct the uncorrelated reference sample, 
taking momentum vectors of protons corresponding to different real events. A real 
event is determined by the momentum vectors of two protons registered in coincidence, 
and an uncorrelated event will thus comprise momentum vectors of protons ejected 
from different reactions. Unfortunately, in such a sample of uncorrelated momentum 
vectors, due to the loss of the kinematical bounds, the production of the r\ meson will 
be not reflected on the missing mass spectrum and hence it cannot be used to extract a 
number of mixed-events corresponding to the production of the r\ meson. Therefore, in 
order to determine a background-free correlation function for the pp — > pprj reaction 
we performed the following analysis: First, for each event, the probability ui was 
determined that this event corresponds to the pp — ¥ pprj reaction. The probability Ui , 
that the i th pp — > ppX event with a missing mass rrii, and a relative momentum of ki 
corresponds to a pp — > pprj reaction was estimated according to the formula: 

w . = Nr,(mi,ki) 

N v (mi,ki) + Np ions (mi,ki)' 
where N v stands for the number of the pp pprj reactions and N p i ons is the number 
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of events corresponding to the multi-pion production with the invariant mass equal to 
m,. The values of N^ijn^k) and N pions (rn, k) were extracted from the missing mass 
distributions produced separately for each of the studied intervals of k. An example of 
a missing mass spectrum with pictorial definitions of N p i ons and N v is presented in the 
middle panel of Figure [2] Now having introduced the weights we can calculate also the 
value of Y*(k) separately for the ppi] and pp(pions) final states. We can achieve this 
by sorting an uncorrelated sample according to the k values similarly as in the case 
of the correlated events and next for each sub-sample we construct background free 
Y*(k) distributions as a sum of the probabilities that both protons in an uncorrelated 
event originate from the reaction where the r/ meson was created. Specifically, if in a 
given uncorrelated event denoted by I, one momentum is taken from a real event say 
11 and the second momentum from another real event 12, then the probability that 
both correspond to reactions where an 77 was created equals to luh ■ w; 2 , and hence 
the uncorrelated yield Y*(k) may be constructed from the sum J2i w zi ' W Z2, where I 
enumerates events in the uncorrelated sub-sample selected for a momentum range k. 
Analogously one can calculate Y* p , iong s assigning to the event a weight equal to 

1-U1. 



2.2. Acceptance corrections 

As the next necessary step in the data evaluation we corrected the determined yields 
to account for the finite geometrical acceptance and detection efficiency of the COSY- 
11 detectors [22]. Hereafter for simplicity by "acceptance" both the geometrical 
acceptance and detection efficiency will be denoted. The acceptance was calculated 
as a function of the proton momentum in the proton-proton rest frame k. It was 
obtained using Monte-Carlo simulations and evaluated according to formula: 

where N acc (k) and N gen (k) denote the number of accepted and generated events, 
respectively. Simulations were performed based on the GEANT-3 packages [23] 
including the realistic geometry of the detectors and a precise map of the field of 
the dipole magnet. The momentum and spatial beam spreads, multiple scattering, 
proton-proton final state interaction [24] and other known physical and instrumental 
effects were taken into account |171fT8] , Knowing the acceptance [22] it would be 
straightforward to correct the nominator of equation [T] however the correction of the 
uncorrelated yield Y* (k) is not trivial since the momenta of protons in the uncorrelated 
event originate from two idependent real events which in general could correspond to 
different values of the detector acceptance. 

Therefore, in order to derive a correlation function corrected for the acceptance, we 
have created a sample of data that would have been measured with an ideal detector. 
For this aim each experimental pp — > ppi] event was multiplied by 1/A(k). This means 
that a given reconstructed pp — > pprj event with momentum of k was added to the 
experimental data sample 1/A{k) times. 

Based on this corrected data sample we calculated the two-proton correlation function 
according to equation [T] In order to avoid mixing between the same events, a "mixing 
step" in the calculations was set to a value bigger than the inverse of the lowest 
acceptance value. The random repetition of identical combinations was also omitted 
by increasing correspondingly the "mixing step". In particular, a \ th event, from 
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the acceptance corrected data sample, was "mixed" with a (\+n) th event, where 
n > max(l/A(k)). If the (l+l) th event was the same as I th , then this was mixed 
with a (l+l+2n) event, etc. 



2.3. Results and conclusions 

The two-proton background-free correlation functions for the pp — > pprj and pp — > 
pp + pions reactions corrected for the acceptances are presented in Figure [3] (full 
squares). As mentioned in the introduction the shape of the obtained correlation 
function reflects not only the space-time characteristics of the interaction volume but 
it may also be strongly modified by the conservation of energy and momentum and by 
the final state interaction among the ejectiles. In order to estimate the influence on the 
shape induced by the kinematical bounds we have constructed the correlation functions 
for both, the pp — » pprj and pp — > pp+pions reaction assuming a point-like source using 
a Monte-Carlo simulation. The results of the simulations are presented in Figure [3] 
(open squares) and it is apparent that they differ significantly from the experimental 
correlation function. At the origin the Coulomb repulsion between protons brings the 
experimental correlation function down to zero. 




k [ MeV/c ] k [ MeV/c ] 

Figure 3. Acceptance corrected experimental proton-proton correlation 
functions for the production of the r) meson (left panel) and multi-pions with 
the invariant mass equal to the mass of the r) meson (right panel). Full squares 
denote the experimental data and open squares represent results of simulations 
determined for a point-like source assuming the homogeneous population of the 
phase space for the reaction products. 



Next, in order to extract from the experimental data the shape of the correlation 
function free from the influence of energy and momentum conservation we constructed 
a double ratio: 

where Const denotes the normalization constant, and indices 'exp' and 'MC refer 
to the experimental and simulated samples, respectively. The determined double 
ratios are given in Table [1] and presented in Figure 0J The double ratio is a well 
established measure of correlations used e.g. by the ALEPH, OPAL and DELPHI 
collaborations for studying the Bose-Einstein or Fermi-Dirac correlations e.g. in the 
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Figure 4. The two-proton acceptance corrected correlation functions normalized 
to the corresponding simulated correlation function for a point-like source. Results 
for the pp — > ppr) (full squares) and pp — > pp + pions (open squares) are compared 
to the two-proton correlation function determined from heavy ion collisions 
(triangles) [25]. 

decays of Z boson 126-23] or W-pairs [53]. In FigureSJa significant discrepancy between 
the two-proton correlation functions determined from inclusive heavy ion reactions 
(triangles) i25j and from exclusive proton-proton measurements (full squares) is clearly 
visible. The data from the kinematically exclusive measurement do not reveal a peak 
structure at 20 MeV/c. Although a similar enhancement is also seen in the invariant 
mass of the proton-proton distributions 20 it disappears in the correlation function 
for the pp — > ppr] and pp — > pp + pions reactions. 

At present it is not possible to draw a solid quantitative conclusion about the size 
of the system since e.g. in the case of the pp — > pprj reaction it would require to solve 
a three-body problem where pp and pi] 30 interactions are not negligible and both 
contribute significantly to the proton-proton correlation. However, based on semi- 
quantitative predictions [21] one can estimate that the system must be unexpectedly 
large with a radius in the order of 4 fm. This makes the result interesting in context 
of the predicted quasi-bound nNN state [32] and in view of the hypothesis [14] that 
at threshold for the pp ~ > pprj reaction the proton-proton pair may be emitted from a 
large Borromean like object whose radius is about 4 fm. 
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Table 1. The double ratios determined for the production of the r\ meson and 
multi-pions with the invariant mass equal to the mass of the r\ meson. In the first 
column of the table k denotes the center value of a 2.5 MeV/c interval. 
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